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^ 1. Introduction 



Abstract. One of the major problems concerning cooling flows 
is the nature of the mechanism powering the emission lines of 
the optical filaments seen in the inner regions of cooling flows. 
In this work we investigate the plausibility of Alfven heating 
(AH) as a heating/excitation mechanism of optical filaments in 
cooling fiows. We use a time-dependent hydrodynamical code 
to follow the evolution of cooling condensations arising from 
the 10^ K cooling fiow. The filaments contain magnetic fields 
and AH is at work at several degrees of efficiency (including 
no AH at all) . We consider two damping mechanisms of Alfven 
waves: nonlinear and turbulent. We calculate the optical line 
emission associated with the filaments and compare our results 
to the observations. We find that AH can be an important 
ionizing and heating source for class I filaments in the line-ratio 
scheme of Heckman et al. In addition, AH can be an important 
contributor to [OI]A6300 emission even for the more luminous 
class II systems. 

Key words: magnetic fields - turbulence - galaxies: clusters 
- cooling flows - intergalactic medium - X-ray: galaxies 



The short cooling times (less than a Hubble time) of the cen- 
tral intracluster medium (ICM) of X-ray emitting clusters of 
galaxies points to gas sinking towards the center of the cluster 
in a "cooling flow" (for reviews, see Fabian et al. 1984,1991 
and Fabian 1994). A sizeable proportion of the cooling flows 
also contains extended (from < 1 kpc to tens of kpc) optical 
emission line filaments around the central dominant galaxy of 
the cluster, which is always at the center of the flow (Heckman 
1981; Cowie et al. 1983; Hu et al 1985; Johnstone et al. 1987; 
Heckman et al. 1989; Baum 1992). These extended optical fil- 
aments systems are seen only in cooling flows and are thought 
to be a phase of the evolution of thermal instabilities arising in 
the cooling flow. The line emission flux ratio [NII]A6583/Hq: 
allowed Heckman et al. (1989) to divide the filaments into 
two distinct classes, class I with an average ([NII]/Hq) — 2.0, 
and class II with {[Nil] /Ha) = 0.9. Also, the class II filament 
systems tend to have higher Ha luminosity and to belong to 
cooling flows with greater mass accretion rates. More recently, 
however, the discovery of fllaments in cooling flows whose line 
ratios are intermediate between those typical of class I and II 



Send offprint requests to: A. C. S. Friaga 



systems (Crawford & Fabian 1992; Allen et al. 1992; Craw- 
ford et al. 1995) called for a continuous distribution instead 
of a clear split into two classes. In addition, there are some 
Ha-luminous class I (i.e. with high [NII]/Ha) systems (e.g., 
A1068, A2146, RX J0439.0-I-0520) departing from the trend of 
class I systems having low Ha luminosities. One of the major 
problems associated with the optical fllaments is their high lu- 
minosities (typically Ha luminosities are 10*^ to 10*^ ergs s"^). 
Quiescently cooling gas that recombines only once at the X- 
ray determined mass accretion rate would not be detectable as 
emission line nebulae. The fact that the central regions of cool- 
ing flows often do show optical emission requires some source 
of energy to repeatedly reionize the gas. Several models, con- 
sidering different mechanisms for ionization and heating of the 
gas, - shocks (David et al. 1988, David & Bregman 1989), 
thermal conduction (Bohringer & Fabian 1989), and photoion- 
ization by soft X-rays and EUV produced in the cooling gas 
(Voit & Donahue 1990; Donahue & Voit 1991; Voit et al. 1994) 
-, have met difficulties in explaining the observed line ratios 
and luminosities. It is possible that some combination of pho- 
toionization and shock models could explain the observations. 
One example of these hybrid models involves the combination 
of self-absorbed irradiating mixing layers of cold clouds embed- 
ded in the hot cooling gas and emission from shocks generated 
in collisions between these clouds (Crawford & Fabian 1992). 
Jafelice and Friaga (1996) have investigated magnetic recon- 
nection as a heating mechanism and concluded that magnetic 
reconnection cannot be invoked as the sole mechanism pow- 
ering the optical line emission in cooling flows. However, it 
could be an important ingredient to explain the emission com- 
ing from low ionization lines. The models with high magnetic 
reconnection efficiencies exhibit a strong [OI]A6300 emission, 
and, therefore, magnetic reconnection could be an additional 
component together with other mechanisms producing most of 
the emission in Ha and higher ionization lines (as [NIIJA6583). 
The existence of magnetic flelds in the ICM has suggested mag- 
netic reconnection as a heating mechanism. In this paper we 
investigate another heating mechanism expected to be at work 
in the presence of magnetic fields: Alfven heating (AH). Here 
we explore the importance of AH not only as a dominant en- 
ergy source for optical filaments in cooling fiows but also as an 
ingredient contributing to the emission of some optical lines. 
The Alfven heating comes from the dissipation (damping) of 
Alfven waves. These waves are easily generated in many cos- 
mic plasmas, but they possess no linear damping mechanism 
since they are not compressive. Nonlinear damping of these 
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waves occurs when one Alfven wave decays into another plus 
a slow magnctosoriic wave, or two Alfven waves combine into 
one fast magnetosonic wave; the resulting magnetosonic waves 
can then be dissipated. The dissipation of these waves may con- 
tribute to heat the solar corona, to influence the interplanetary 
wave spectrum, and to determine the wave spectrum available 
to scatter cosmic rays. As they are not compressive, the non- 
linear decay occurs due to second order effects in the wave 
amplitude. Interacting Alfven waves are compressive and thus 
generate a new compressive wave (see Chin and Wentzel 1972). 
The resulting heating — at least in terms of nonlinear, reso- 
nance surface and turbulent damping mechanisms — has been 
applied to study the thermal stability in the broad line region 
of quasars and in the basis of the wind in hot stars, in a suc- 
cessful way (Gongalves et al. 1993a, 1996a,b). Alfven waves 
are created by perturbed magnetic fields. In view of this, two 
conditions that are fulfilled by the core of a cooling fiow make 
it a favorable site for the generation of Alfven waves: it con- 
tains magnetic fields, and it is turbulent. The lack of any strong 
hard X-ray emission due to inverse Compton from the relativis- 
tic electrons in a few clusters (generally not cooling fiows) that 
have halo radio sources imposes a lower limit for the magnetic 
field of 0.1 /iG over scales of up to 500 kpc (Rcphacli & Grubcr 
1988). At the smaller (~ 100 kpc) scales of cooling flows, mag- 
netic fields of typically 1 — 3 /iG are derived corresponding to 
a magnetic pressure of about 1% of the thermal pressure. At 
the center of a cooling flow still more intense magnetic flelds 
must be present since the fleld is amplified by compression as 
the gas fiows inward (Soker & Sarazin 1990). As a matter of 
fact, Faraday rotation and depolarization of the radio emis- 
sion of extended radio sources (Ge & Owen 1993; Taylor et 
al. 1994) residing in cooling flows have revealed that the mag- 
netic field increases inward and that it can reach 20-100 fj,G 
depending on the degree of ordering of the field. In this way, 
at the very center of the fiow, the magnetic pressure can su- 
persede the thermal pressure. One kinematic signature of the 
filament systems is the lack of organized velocity patterns (e.g., 
rotation or shear or infall) (Heckman et al. 1989; Baum 1992). 
The filaments typically have very small rotational velocities 
and are turbulently, and not rotationally, supported. A small 
number of filament systems do show apparent rotation with 
rotational velocities of ~ 300 km s~^ within a few kpc. How- 
ever, these same systems show disordered velocity patterns at 
larger scales (Hu et al. 1985). Lines are broad throughout the 
filamentary region, but line widths decrease from 500-100 km 
s~^ in or near the galactic nucleus to 100-300 km s~^ at the 
largest radii (typically 5-15 kpc). This suggests that the inner 
ICM is highly turbulent, with the amplitude of random veloci- 
ties increasing inward. In fact, the core of a cluster is expected 
to be turbulent not only due to the frequent crossing of galaxies 
through it but also due to relics of merging of subclustcrs. In 
addition, the absence of rotation in the center of cooling fiows 
requires rotational breaking of the gas flowing inward, and tur- 
bulent viscosity is the most likely transport process of angular 
momentum in a cooling flow (Nulsen et al. 1984). The possi- 
bility of using the turbulent energy to heat the plasma in the 
core of cooling fiows was also considered by Loewenstein and 
Fabian (1990), in order to explain the kinematics of the cooling 
flow clouds. They assume a magnetic viscous heating process 
(' plasma slip') where the magnetic fleld that passes from the 
clouds to the hot gas can be forced to oscillate by the noise 
and so cause the ionized particles within the cloud to oscillate 



and collide with neutral particles. This process can efficiently 
transport Alfven wave energy to length scales of ~ 10^^ cm. 



2. Alfven Heating in the Intracluster Medium 

The emission line filaments have line widths of several hundred 
to -1000 km s"^ (Heckman et al. 1989; Baum 1992; Crawford 
& Fabian 1992; AUen et al. 1992; Crawford et al. 1995). It is 
possible that these velocities of the filaments are turbulent and 
that the hot component of the ICM is itself turbulent. Anyway, 
there is plenty of kinetic energy to be tapped via Alfven waves 
into an energy source for the ionization of the filaments. Several 
damping mechanisms can be responsible for Alfven heating. In 
this work we consider the nonlinear damping and the turbulent 
damping for Alfven waves. The damping mechanisms above 
have been investigated before in many astrophysical objects: 
late-type stars, protostcUar and solar winds (Jatenco-Pereira 
& Opher 1989a,b,c); galactic and extragalactic jets (Opher & 
Pereira 1986; Gongalves et al. 1993b); early-type stars (dos 
Santos et al. 1993); broad line regions of quasars (Gongalves et 
al. 1993a, 1996a) and others. 



2.1. Nonlinear damping 

An Alfven wave is likely to dissipate because of its nonlin- 
ear interaction with either the non-uniform ambient field or 
another Alfven wave (Wentzel 1974). The nonlinear interac- 
tion of maguetohydrodynamic waves has been treated in de- 
tafl by Kaburaki & Uchida (1971), Chin & Wentzel (1972), 
and Uchida & Kaburaki (1974). When the magnetic fleld is 
weak, i.e. when the Alfven velocity va is lower than the sound 
velocity Cs, it is found that two Alfven waves travelling in op- 
posite directions along a magnetic fleld line can couple nonlin- 
early to give an acoustic wave, which in turn dissipates rela- 
tively quickly. In regions of strong {va > Cs) magnetic field, 
one Alfven wave can decay into another Alfven wave and a 
sound wave travelling in the opposite direction. The result- 
ing Alfven wave has a frequency smaller than the original one 
and it can, in turn, decay into another lower-frequency Alfven 
wave plus an acoustic wave. The cascade continues until all 
the Alfvenic energy has been converted to acoustic waves that 
dissipate rapidly. The latter possibility is considered in this pa- 
per. In our models, the onset of nonlinear heating occurs for 
(3 = Ps/Pgas = iB^/8n)/nkBT > f3o„ > 1, where Pon is the 
initial value of f3 when AH becomes important (/3 > 1 corre- 
sponds to va/cs > \/G/5 for an ideal 7 = 5/3 gas). Note that, 
following the convention widely used in cooling fiow studies, we 
define /3 — Pb / Pgas , whereas in plasma physics, l3pi = Pgas/Ps 
denotes the "beta parameter" of the plasma. According to La- 
gage & Cesarsky (1983), the nonlinear damping rate is 

^~^-4V2^'"W''sV^ ' 

where ^ = 5 — 10, Cs is the sound velocity, va the Alfven veloc- 
ity, Lu a characteristic Alfven frequency and p{Sv'^) / (B^ /8tt) 
the ratio of the energy density of Alfven waves to that of the 
magnetic fleld. As we are interested in AH coming from the 
nonlinear damping of Alfven waves, we heuristically derive the 
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dependences of the heating rate Hnl (erg cm ^ s ^) with den- 
sity and temperature. Thus, the nonhnear Alfvenic heating is: 



VA 



where 3>u, is the wave flux, Tnl is given by (1), and va = 
Bj^f^np is the Alfven velocity. As our perturbation collapses, 
the component of the magnetic field perpendicular to the di- 
rection of compression is amplified. For our plane parallel ge- 
ometry, B oc (X p. where A is the cross-sectional area per- 
pendicular to the magnetic field. As a consequence, va oc p^^"^ ■ 
Let the dependence of $to on p be given as oc p'^ . Taking 
= p{5v'^)vA, we have p{&v ) oc p"'^'"^. Since Cs oc T^/^, for 
w independent of the density, we obtain: 



Hnl oc p 



2C-7/2 rpl/2 



(3) 



For a plane-parallel geometry, oc ^ ^ oc p, implying c = 1, 
and equation (3) becomes: 



Hnl oc p-^/^ t"^ 



(4) 



2.2. Turbulent Damping 

Turbulence theory indicates that the spectrum of turbulent 
Alfven waves has a power law dependence in phase space 

P{k) oc , 

where k is the Alfven wavenumber, p is the spectral index, and 



P{k)4'n-k^dk 



is the energy density of turbulent Alfven waves. The value of 
the spectral index is determined by the type of the turbulence 
(De Groot & Katz 1973; Kainer et al. 1972; Eilek & Henrik- 
sen 1984). Tu et al. (1984) noted that the nonlinear interac- 
tion between outward and inward propagating waves results 
in an energy cascading process. A small part of the energy 
of lower frequency fluctuations cascades to higher frequency 
fluctuations until it reaches a frequency high enough for some 
dissipation process to occur, the most probable candidate be- 
ing proton-cyclotron damping. As noted by Hollweg (1987), 
the observed magnetic power spectra, Pb, in the solar wind 
exhibits a well- developed inertial subrange resembling Kol- 
mogorov turbulence, i.e. Pb oc k~^^^ , over more than 3 orders 
of magnitude of k (Matthaeus & Goldstein 1982). This indi- 
cates a cascade to small scales. The turbulent cascade transfers 
wave energy from large to short scales, and, in the end, the dis- 
sipation of energy through turbulence is governed by the small 
scale linear absorption mechanism. Exploiting the similarity of 
Pb oc k~^^^ and Kolmogorov turbulence in ordinary fluids, the 
volumetric heating rate associated with the cascade process 
can be written as 



Ht = p{6vy^/Lc 



(5) 



where p is the mass density, (<5w^) is the velocity variance asso- 
ciated with the field, and Lcorr is a measure of the transverse 
correlation length (Hollweg 1986, 1987). As in the previous sub- 
section, we can derive the dependence of the turbulent heating 



rate on p. For a plane parallel geometry, {Sv'^) oc p'^ (as- 
suming ^u> oc p"), and adopting Lcorr oc p'', we obtain: 



(2) Ht oc p(3/2)c-6-5/4 



(6) 



For a plane parallel geometry, c = 1 ($to oc p), and, if we as- 
sume that Lcorr oc B'^^^ (Hollweg 1986, 1987), then 6 = -1/2 
(since B cc p). The resulting turbulent Alfven wave heating is 



Ht oc p 



3/4 



(7) 



3. Time-dependent Models for Cooling Condensations 

We have investigated the evolution of optical filaments since 
their formation out of the hot phase of the cooling flow and 
calculated their optical signature within the scenario outlined 
in Section 2. The evolution of the cooling filaments is obtained 

by solving the liydrodynamical equations of mass, momentum 
and energy conservation (see Friaga 1986,1993; Friaga & Ter- 
levich 1994,1996). The condensations are assumed to be self- 
gravitating. Since there is no ionization equilibrium for tem- 
peratures lower than 10® K, the ionization state of the gas at 
T < 10^ K is obtained by solving the time-dependent ioniza- 
tion equations, for all ionic species of H, He, C, N, O, Ne, 
Mg, Si, S, Ar and Fe. We adopt an isochoric non-equilibrium 
cooling function for temperatures lower than 10^ K, since the 
recombination time of important ions is longer than the cool- 
ing time at these temperatures. We have chosen an isochoric 
cooling function because the evolution of the gas is nearly iso- 
choric as soon as the gas cools below 10® K, first due to the 
rapid cooling for temperatures around the peak (T ~ 10^ K) 
of the cooling function, and then due the magnetic support 
against further compression of the gas at lower temperatures. 
The cooling function and the coefficients of coUisional ion- 
ization, recombination and charge exchange of the ionization 
equations are all calculated with the atomic database of the 
photoionization code AANGABA (Gruenwald & Viegas 1992). 
The adopted abundances are half-solar, as derived from X-ray 
spectroscopy of the ICM. (Solar abundances are from Grevesse 
& Anders 1989.) The initial density perturbations are char- 
acterized by an amplitude A and a length scale L, following 
5p/p = Asin{x)/x (here x = 2i:r/L). We have also assumed 
a plane-parallel geometry and that the perturbations are iso- 
baric and nonlinear {A = 1). The geometry is justified by the 
fact that the line emission is filamentary, which suggests that 
the perturbations are sheetlike rather than spherical. We start 
to follow the evolution of the perturbations from the nonlinear 
stage in view of the uncertainties about processes suppressing 
the growth of thermal instabilities in cooling flows (see Sec- 
tion 5). An unperturbed uh = 0.1 cm" and T = 10'^ K IS 
assumed for the filament, appropriate for the inner regions of 
a cooling fiow, where the filaments are more often seen. L was 
fixed at 1 kpc for all the models. This length scale is suggested, 
for instance, by the spatial fluctuations in the velocity of fila- 
ments resolved in nearby cooling flows (Heckman et al. 1989). 
We fixed /3 — Ps/Pgas = 0.1 for the unperturbed medium, a 
representative value for the range of /3 = 0.01 — 1, expected in 
the central 10 kpc of cooling fiows. 

Here we investigate several representative models with mag- 
netic fields, one with no AH and the others with AH for two 
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Table 1. Properties of models. 



Model 


Heating 


c 




tcol 


tmax 


ii/2 


L H a.max 


L H a,ave 










(10'' yr) 


(10" yr) 


(10" yr) 


(10"° erg s"^) 


(10*° erg s 


A 


none 


0. 




6.00 


0.345 


4.51 


2.65 


1.08 


B 


turbulent 


1. 


1. 


7.36 


4.56 


4.10 


39.8 


27.1 


C 


turbulent 


0.3 


1. 


6.24 


0.686 


4.51 


4.30 


2.81 


D 


nonlinear 


0.3 


3. 


6.16 


0.467 


4.41 


3.19 


1.44 


E 


nonlinear 


1. 


10. 


6.00 


3.67 


4.51 


20.9 


7.70 


F 


nonlinear 


0.3 


10. 


6.00 


0.447 


4.51 


3.20 


1.40 



damping mechanisms (nonlinear and turbulent) and at differ- 
ent efficiencies (see Table 1). Model A allows for the presence 

of a magnetic field only through a magnetic pressure term 
(Pb = /8it) in the equation of motion. It gives the bench- 
mark to evaluate the effects of AH. In models B to F, the 
energy equation include AH as an additional heating term, 
which is turned on only for /3 > f3on, where /3on is the value of 
/? when AH becomes important. Wo considered nonlinear and 
turbulent Alfven wave heating. From equations (4) and (7), 
the heating term has the form 

Hnl = CHoHnoy'^^T/ny^^ , (8) 
for nonlinear heating, and 

Ht = C Ho{n/nof^^ , (9) 

for turbulent heating, where is the efficiency of Alfven heat- 
ing. The choice of the normalizations is the following: 1) 
To = 10^ K (the temperature at which optical emission be- 
gins to be important); 2) no = 3.267 x 0.1 cm~^, where 3.267 
is the compression factor from the unperturbed uh = 0.1 cm~^ 
and T = 10^ K initial state, to T = To, under isobaric con- 
ditions; and 3) Hi, = (3/2) nksT/tcoi = 2.55 x 10"^" erg s"^ 
cm""*, where tcoi is the collapse time (defined below) for the 
filament with no heating. 

4. Results 

In order to assess the importance of the optical emission, a 
number of quantities characterizing the optical phcise of our 
models are presented in Table 1: the collapse time tcoi is the 
time elapsed since the beginning of the calculations until the 
beginning of the optical phase; LHa,max is the maximum lu- 
minosity in Ha; tmax is the time when LHa.max is attained; 
ti/2 is the time when half the mass of the perturbation has 
cooled below T = 5 x 10® K; [tmax and ti/2 are counted from 
the beginning of the optical phase); and LHa,ave is the average 
luminosity since the beginning of the optical phase until ii/2. 
The luminosity Lho is normalized to M = 100 M© yr~^, with 
M taken as equal to the ratio of the mass of the condensation 
to its collapse time. In the beginning of the growth of the per- 
turbation, the magnetic pressure is unimportant in comparison 
to the thermal pressure. However, as the gas within the per- 
turbation cools, under conditions of a frozen-in field, the com- 
ponent of the magnetic field perpendicular to the direction of 
compression increases (S oc p in plane-parallel geometry), and 
eventually (3 becomes > 1. At this stage, the magnetic pressure 



Fig. 1. Evolution of Lua, the normalized Ha luminosity (for its 
definition see the text), for the models A-C (solid lines), D (dotted 
line) and E-F (dashed lines); the curves are displaced upwards for 
increasing values of C,. 

halts further compression, while the temperature keeps drop- 
ping (if no heating mechanism is efficient enough to prevent 
the gas from cooling). When the temperature has decreased 
enough, the gas begins to emit optical lines. We consider the 
beginning of the phase of optical line emission when the tem- 
perature of the innermost cell (i.e. the first cell to cool) falls 
below 5 X 10® K. Before this temperature has been reached, 
the gas is in its X-ray phase, during which most of the cool- 
ing occurs via EUV and X-ray lines. The collapse time is the 
smallest (6 x 10® yr) for models with no heating (model A) or 
for those in which the heating is turned on later in the evo- 
lution of the perturbation (models E and F). As the heating 
becomes more important during the X-ray phase, the collapse 
is delayed (i.e. tcoi is increased). Model B (turbulent heating 
and = 1) is the slowest to collapse. A model with nonlin- 
ear heating, /3on = 3, and C = Ij has been run, but its mass 
that reaches T < 5 x 10® K is at most 35% of the total mass 
of the perturbation. This result suggests that nonlinear Alfven 
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heating can constitute a mechanism for suppressing the growth 
of thermal instabihtios in cooUng fiows. The nonhncar and the 
turbulent AH's represent extreme opposite dependencies of AR 
on density. Turbulent AH {Ht oc n^^*) deposits more energy 
at the center of the filament, where the density is higher. On 
the other hand, nonlinear AH {Hnl oc n~^^'^T^^^) deposits 
more energy per unit volume in the outer regions of the fila- 
ments, whore the density is lower. However, due to magnetic 
support, the central density of the cooled filament (T « 10^ 
K) is only « 3.3 times the initial unperturbed density, and 
the central decrease of the AH rate due to density variation 
amounts to only a factor 6. In the end, the fact that the en- 
ergy is deposited in a much cooler gas at the center than on the 
outside of the filament (T « 10^ K) makes the AH somewhat 
more important in the colder, central gas than in the outer skin 
of the filament, as we can see from the timescale for heating 
tn = {3/2)nkBT/HNL oc n^^^T^/^, which is 1.6 times smaller 
at the center than on the outside of the filament. The inverse 
dependency of Hnl on the density explains the suppression of 
the growth of the perturbation in the model with non-linear 
heating, C = 1 and 0on = 3. It should be noted that in our 
calculations, we assume that the filaments are optically thin to 
Alfvcn waves, so that energy is deposited by AH throughout 
the filament. The evolution of Lho is shown in Figure 1. and 
the evolution of the line ratios can be followed in the diagram 
[NII]A6583/Ha versus [OI]A6300/Hq (Fig. 2). For a model to 
be successful, it should not only reproduce the position of the 
filaments in the line ratio diagrams but also account for their 
luminosity. The models B (turbulent heating) and E (nonlin- 
ear heating) are the more energetically promising. The energy 
budget requirements of class I filaments are satisfied (see Ta- 
ble 1 and Fig. 1) since the class I objects arc the less luminous 
systems (an average (Lhc) = 4 x lO*" erg s^^). Our models 
oven fulfill the luminosity requirements for some of the less lu- 
minous class II filaments, with an average {Lho) = 3 x 10'*^ 
erg s~^ (Heckman et al. 1989). The luminosity of our mod- 
els, however, fall short of the most luminous class II systems 
(e.g. Perseus), so that we propose the AH as the mechanism 
powering primarily class I filaments. From the [Nil] /Ha ver- 
sus [01] /Ha diagram, the model B (with turbulent heating) 
reproduces the values of [OI]/Ha of the filaments, but some- 
what overproduces [NII]A6584 emission of class I filaments. On 
the other hand, model E (with nonlinear heating) reproduces 
the values of [NII]/Ha, but overproduces [OI]A6300 emission. 
It seems that if some combination of these mechanisms is at 
work, a range of heating efficiencies and chemical abundances 
could account for the loci of class I systems in the [Nil] /Ha 
versus [OI]/Ha diagram. 



5. Conclusions 

From the energetic viewpoint, our models could account for 
typical filament systems with Hmc « 10 — 100 but not for the 
extreme ones (e.g., Perseus) with Hrec ~ 1000 {Hrec is the 
number of recombinations per cooling proton required to yield 
the observed Ha luminosity). Our models are better at repro- 
ducing the luminosities of the weaker class I filaments than the 
more luminous class II systems. AH could be a dominant heat- 
ing mechanism for the class I filaments, and even for the less 



luminous class II filaments. It should be noted that there are 
evidences that matter is dropping out of the cooling fiows at 
radii larger than r ~ 10 kpc, so that only a fraction of the over- 
all cooling flow rate (over typical radii of ~ 100 kpc) reaches the 
region where the filaments are seen. A fraction (f) w 0.1 reach- 
ing r ~ 10 kpc is suggested from a M oc r variation (Thomas, 
Fabian & Nulscu 1987). This value is derived from a simple 
model in which the flow is stationary and the X-ray emission 
comes from a homogeneous gas. Taking into account nonsteady 
flow (the flow time from 100 to 10 kpc is comparable to the age 
of the system — for instance, the values of n and T adopted in 
the present calculations imply a flow velocity of 5 — 50 km s~^ 
at 10 kpc for M = 10 - 100 Mq yr"^) and the contribution of 
nonlinear blobs to the X-ray emission makes the reduction of 
M less drastic than that given by the M oc r model (Meiksin 
1990; Friaga 1993). These effects, however, do not prevent mat- 
ter from being removed from the flow, and the resulting (f> is 
~ 0.2 — 0.3. Therefore, mass removal from the cooling flow lim- 
its the energetic input of AH, and the luminosities given in the 
previous section should be reduced by a factor 0. Even if AH 
cannot be invoked as the sole mechanism powering the opti- 
cal line emission of class II fllaments, it could be, however, an 
important ingredient to explain the omission coming from low 
ionization lines in these systems. The models with nonlinear 
heating and high AH efficiencies exhibit a strong [OI]A6300 
emission, and, therefore, AH could be an additional compo- 
nent together with other mechanisms for producing most of 
the emission in Ha and higher ionization lines (as [NIIJA6583). 
The combination of two distinct heating mechanisms, each one 
being responsible for certain emission lines and/or explaining 
one class of filaments, has allowed to build hybrid models for 
optical filaments, which have been particularly succesful at re- 
producing the observations of optical filaments. Shocks models 
tend to reproduce the locus in the line-ratio diagram of class 
II filaments, although the energetic requirements arc stringent. 
A hydrid model could, therefore, combine shocks and an ad- 
ditional heating source which explains the higher ionization 
class I filaments. In this way, Crawford & Fabian (1992) in- 
voke shocks between cold clouds to explain the class II fila- 
ments (and produce their [01] emission) and mixing layers to 
account for the class I filaments with lower Ha luminosities 
and higher levels of ionization (higher [NII]A6583/Ha ratios). 
Other additional heating sources are magnetic reconnection 
(Jafelice & Friaca 1996), and Alfvcn heating (this work). Jafe- 
lice and Friaga (1996) suggest a hydrid model involving mag- 
netic reconnection (with an efficient producion of [OI] emis- 
sion), and self-absorbed mixing layers with high flux and low 
temperature for the incident radiation (producing most of the 
Ha and [Nil] emission). Whether AH is at work together with 
shocks or with mixing layers, it efficiently supplies [01] emis- 
sion. For a representative Lna = 1 x W^^(j) erg s^^ of model 
E, one obtains a [OI]A6300 luminosity « 2 x W^^cj) erg s~^ (see 
Fig. 2). Given the average ratio [OI]A6300/Ha = 0.25 of the 
optical filaments, AH could account for the [OI]A6300 emis- 
sion of systems with Lna up to « 8 x W^^cj> erg s^^. The AH 
keeps the gas relatively warm, giving rise to strong emission of 
[OI]A6300. We note that the energetic requirements concerning 
the [OI] emission of the more luminous systems, like Perseus, 
with Lhci = 2.6 X lO''^ erg s~^ (considering a Ha luminosity 
of 4.7 X lO''^ erg s~^ (Heckman et al. 1989) and a cooling flow 
rate of 183 Mg yr^^ (Fabian 1994)), arc not satisfied by our 
model even in the most favorable case. For these extreme sys- 
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tems, an additional heating mechanism is required. Whereas 
there is undoubtedly a connection between cooling flows and 
optical filaments, this connection is troubled by some facts, 
and our model can shed light in this issue. In the first place, 
while there is a strong correlation between the Ha luminosity 
and the mass accretion rate (at a 99.94% confidence level ac- 
cording to Hcckman ct al. 1989), the relation cannot be direct 
because of the large (more than two orders of magnitude) scat- 
ter in the relation Lhoi — M and the fact that some massive 
cooling flows (e.g., A2029) show no detectable line emission. 
The existence of a cooling flow, therefore, is a necessary but 
not sufficient condition for detectable line omission. Secondly, 
except for a few filament systems (e.g., Perseus and A1795) 
extending over several tens of kpc, the optical line emission is 
confined to the inner 10 kpc of cooling flows, whereas the X- 
ray images of clusters reveal that the cooling gas is distributed 
throughout the cooling radius (~ 100 kpc). The high spatial 
concentration of optical filaments towards the center could be 
explained by a greater magnetic field strength in the central 
regions (Soker & Sarazin 1990) and by a higher level of turbu- 
lence in these regions (Loewenstein & Fabian 1990; Lowewen- 
stein 1990; Begelman & Fabian 1990). As a consequence, an in- 
tense Alfven wave generation is boosted only in the inner cool- 
ing flow. On the other hand, the level of turbulence could be 
the "second parameter" besides mass accretion rate, regulating 
the luminosity of optical filament systems. Turbulence, caused, 
for instance, by stirring of the central ICM by a galaxy moving 
through the cluster core, or by interactions between the ICM 
and the relativistic plasma of an extended radio source, could 
trigger, via generation and dissipation of Alfven waves, optical 
emission lines. In this scenario, the generation of Alfven consti- 
tutes a step in the tapping of the turbulent energy of the ICM 
into line luminosity. Additional support for the idea that tur- 
bulence can be crucial for producing the line luminosity comes 
from the fact that the kinetic energy flow, determined from the 
observed velocity dispersion of the filaments, correlates much 
tighter (a scatter of one order of magnitude) with the Ha lumi- 
nosity (Heckman et al. 1989) than M does. This indicates that 
it is the combination of mass flow and high levels of turbulence 
that is relevant for the line luminosity, rather than a large mass 
accretion rate alone. The present calculations are intended to 
apply to the inner ~ 10 kpc of the cooling flows, where the 
optical filaments are most commonly seen. The turbulent ve- 
locities of the cooling flow increase toward the center of the 
cluster, as inferred from the increase of the line widths near 
the center (Heckman et al. 1989). Since the turbulent energy 
of the cooling flow is the energy source of the Alfven waves, it 
is only in the inner ~ 10 kpc that the density of turbulent en- 
ergy is high enough to make AH efficient. In this region, blobs 
condensating out of the cooling flow are kept warm enough 
by AH to produce strong line emission. Blobs cooling down 
in the outer regions are subject to low levels of AH and cool 
quiescently without noticeable optical emission (similarly to 
model A). In this way, the AH mechanism is consistent with 
the fact that line emission does not extend to ~ 100 kpc. A 
further issue raised by the results of our models concerns the 
growth of thermal instabilities in cooling flows. It has been sug- 
gested that a number of processes could inhibit the growth of 
thermal instabilities in the presence of gravitational fields and 
background flow in cooling flows (Malagoli et al. 1987; Balbus 
1988; Balbus & Soker 1989; Tribble 1989; Loewenstein 1989; 
Brinkmann et al. 1990; Hattori & Habe 1990; Yoshida et al. 



1991; MalagoH et al. 1990; Reale et al. 1991). The fact that one 

of our models (with nonlinear heating, f3on ~ 3, and ^ = 1) has 
reached thermal stability shows that also Alfven heating can 
constitute a mechanism for suppressing strong growth of ther- 
mal instabilities in cooling flows. On the other hand, Gongalves 
et al. (1993a; 1996a) have considered the possibility that the 
broad line regions of quasars are formed via thermal instability 
in the presence of AH. They have investigate three damping 
mechanisms of Alfven waves - resonance surface, nonlinear, 
and turbulent -, and found that AH could establish a stable 
two-phase medium in the broad line regions of quasars. The 
thermal stability of the model mentioned above suggests that, 
also in the case of the cooling flow medium, AH could induce 
a stable two-phase equilibrium. 
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